Stable demographic ratios of haploid gametophyte to diploid sporophyte abundance in macroalgal populations

Macroalgal populations often consist of free-living haploid (gametophyte) and diploid (sporophyte) stages. Various ecological studies have been conducted to examine the demographic diversity of haploid-diploid populations with regard to the dominant stage. Here, I relaxed the assumption of classical research that the life history parameters of haploids and diploids are identical and developed a generalized haploid-diploid model that explicitly accounts for population density dependence and asexual reproduction. Analysis of this model yielded an exact solution for the abundance ratio of haploids to diploids in a population in which the ratio is determined by the balance of four demographic forces: sexual reproduction by haploids, sexual reproduction by diploids, asexual reproduction by haploids, and asexual reproduction by diploids. Furthermore, the persistence of a haploid-diploid population and its total biomass are shown to be determined by the basic reproductive number (R0), which is shown to be a function of these four demographic forces. When R0 is greater than one, the haploid-diploid population stably persists, and the ploidy ratio obtained by the analytical solution is realized.


Introduction
Sexual reproduction in eukaryotes leads to the alternation of two nuclear phases, namely the haploid (n) and diploid (2n).In many plants, algae, and fungi both ploidy stages develop into adult free-living organisms and as such, a biphasic life cycle is observed [1][2][3][4][5][6][7].Many macroalgae alternate between a free-living haploid generation (the gametophytes) and the free-living diploid generation (the sporophytes).
The life cycle diversity has been interesting in the context of how the evolution from haploidy to diploidy originally occurred [8,9].For example, the hypotheses that the diploids have an advantage in masking deleterious mutations [9][10][11], that the diploids evolve faster [12,13], and that the haploids have an advantage in poor-nutrient conditions [14][15][16], have been considered.
Population geneticists have focused primarily on life cycle evolution, whereas many phycologists have been interested in the descriptive study of the macroalgal life cycles.Classically, the type of alternation of generations between two distinct free-living stages is classified as heteromorphic (morphologically distinct from one another) or isomorphic (seemingly identical to each other) [17].Species with the heteromorphic life cycle show large phenotypic differences (e.g., macroscopic versus microscopic, erect versus crustose).These species also often show different seasonal appearances and habitats [18].In contrast, the two stages of species with isomorphic life cycle tend to be observed simultaneously and sympatrically [19][20][21][22][23].
Because these studies were based on laboratory observations, there was interest in field research.Especially in species with the isomorphic life cycle, gametophytes and sporophytes (tetrasporophytes in red algae) are often observed simultaneously, and their abundances were investigated.Field studies of the ratio of haploid to diploid abundance have shown both haploid-dominant [20,21,23] and diploid-dominant [19,22] populations exist.Different patterns of haploid or diploid dominance have been observed for different populations of the same species, and even over the same population depending on the special location or season [24][25][26][27].
Theoretical models have been developed to explain the variety of haploid versus diploid dominance patterns [28,[30][31][32][33][34][35][36][37][38].Many studies are based on numerical calculations, and few studies derive simple analytical solutions.In particular, the classical study of Thornber and Gaines [28] presented an analytical solution only for the case where haploids (gametophytes) and diploids (tetorasporophytes) have identical survival and growth rates (i.e., ploidy abundance ratio = ffi ffi ffi 2 p : 1).Although in their classical study, Littler et al. [39] reported the functional similarity in physiological or ecological performances between gametophyte and sporophyte in an isomorphic species, many recent empirical studies have reported various types of functional differences between haploids and diploids.Specifically, chemical composition [40], resistance to physical stress [41], fecundity [22,28], survivorship [29,42,43], resistance to predators [44], resistance to epiphyte infection [45], growth rate [46][47][48], dispersal ability of the reproductive cells [49], and viability of the reproductive cells [50] are reported.Based on these facts, this article analyses a mathematical model that explicitly considers the differences in life history parameters of haploid gametophytes and diploid sporophytes, and aims to express stable ploidy ratios in a haploid-diploid population using a mathematical expression as simple as possible.
In recent years, in addition, the importance of asexual reproduction, which deviates from the normal reproductive process in which gametophytes reproduce sporophytes and sporophytes reproduce gametophytes (see next section), has been pointed out [51][52][53].Unfortunately, theoretical studies that consider asexual reproduction are limited.Among the few studies, [53] analyze the demographic equilibrium state of a population genetic model that considers sporogenesis, and Vieira and Santos [31,35] analyze the population structure yield by a matrix model that considers clonal reproduction (i.e., vegetative growth of new fronds from the same holdfast).In addition, this study explicitly considers density dependence, which many previous studies have omitted for simplicity.Many previous studies (e.g., [31]) define the ploidy ratio by the right eigenvector in the matrix model, where the density dependence is ignored and the population size diverges to infinity.In reality, the population size should remain finite and the stable ploidy ratio of this population is important.Vieira et al. [38] develop the individual-based models to investigate population dynamics that consider the realistic density dependence, but the results are complex and analysis of a simpler model is required.Analyzing the model, I investigate the ploidy ratio in a haploid-diploid population with asexuality and if the ploidy ratio is in agreement with the results of previous studies when the parameters are constrained accordingly.

Life cycle of macroalgae
In this section, I briefly review the life cycle diversity in macroalgae (see also [3,6,7,17,54]), which is the prerequisite knowledge of the models.A diagram of the reproduction process assumed in the model is shown in Fig 1A.

Alternation of nuclear phases and generations
In sexual eukaryotes, states characterized by one set of chromosomes (1n, haploid) alternate with states characterized by two sets of chromosomes (2n, diploid) through meiosis and zygosis.This "alternation of nuclear phases" often results in the alternation of individuals that exhibit different reproductive patterns during their life cycle, which is referred to as the "alternation of generations".In many macroalgae, haploid multicellular individuals and diploid multicellular individuals appear as part of this alternation of generations (haploid-diploid life cycle).Generally, the former is called "gametophytes" because they reproduce by gametes, and the latter is called "sporophytes" because they reproduce by spores (Fig 1B).
I note that, however, the alternation of generations is not an essential phenomenon in macroalgae.For example, individuals belonging to Fucales, a group of brown algae, are diploid.When they mature, gametogenesis and syngamy occur continuously, and the life cycle is completed as the diploid zygote develops into a multicellular organism again (diploid life cycle).In other words, they do not exhibit alternation of generations.

Heteromorphic and isomorphic life cycles
When we focus on the appearance of macroalgal gametophytes and sporophytes, they are often very different in morphology.Extreme examples include species of Laminariales, in which macroscopic sporophytes alternate with microscopic gametophytes, and the species of Scytosiphonales, in which erect gametophytes alternate with crustose sporophytes.This mode of alternation of generations is called "heteromorphic" alternation of generations.On the other hand, (especially immature) gametophytes and sporophytes are often so similar that they are indistinguishable (e.g., Ulva), and such an alternation of generations is called "isomorphic" alternation of generations [17].The life cycle characterized by the heteromorphic (isomorphic) alternation of generations is called the heteromorphic (isomorphic) life cycle.
I note that the difference between gametophytes and sporophytes is sometimes less extreme.In the brown algae Ectocarpus, for example, both gametophytes and sporophytes are small and filamentous, but there is a substantial morphological difference between them.This species is classified as "heteromorphic (near-isomorphic)" [51].The accumulation of empirical studies showing that gametophyte and sporophyte (or tetrasporphyte of red algae) phenotypes differ (morphologically, ecophysiologically, or biomechanically) calls into question the classical classification of the life cycle into heteromorphic and isomorphic categories.

Triphasic life cycle in red macroalgae
When we consider macroalgal life cycle, we generally assume alternation of generations between haploid gametophytes and diploid sporophytes (biphasic life cycle).Among red algae, however, many species from the Florideophyceae are known to exhibit more complex life cycles.Here, syngamy occurs when a haploid sperm without flagella (spermatium) attaches to haploid female gametes (carpogonium) that attaches on the female gametophyte.This fertilized cell develops into diploid multicellular structure, the carposporophyte.First, the carposporophyte produces diploid carpospores by somatic cell division.Only later they are released and these develop into free-living diploid sporophytes, tetrasporophytes.Finally, the tetrasporophyte produces haploid tetraspores by meiosis, and these develop into gametophytes to complete their life cycle.Their life cycle involves three "phases" of haploid gametophyte, diploid carposporophyte, and diploid tetrasporophyte, so this can be called a "triphasic" life cycle [55].
Because the ploidy of carposporphyte (diploid) is different from gametophyte (haploid), the caposporophyte is recognized as a distinguishable phase, but is nutritionally dependent (parasitic) on the female gametophyte [56].Hence, if we consider that carposporophytes have little effect on the ecological properties of the species, the triphasic life cycle of red algae can be approximated as a biphasic life cycle alternating only between gametophyte and tetrasporophyte (Fig 1C).

Asexual reproduction in macroalgae
In classical ecology and evolutionary researches about the macroalgal life cycle, researchers assume obligate sexuality, with gametophytes reproducing sporophytes and sporophytes reproducing gametophytes (e.g., [9,57]).However, many macroalgae exhibit complex reproductive strategies, including asexual reproduction (e.g., [51]).Examples include parthenogenesis, in which gametes that fail to syngamy develop into gametophytes, asexual multiplication of gametophytes by monospores, asexual multiplication of sporophytes by mitospores, and vegetative growth (e.g., the growth of new fronds from the same holdfast or from broken fragments) (Fig 1D ).
For simplicity, this article models asexual reproduction, as processes in which haploid gametophytes reproduce haploid gametophytes and diploid sporophytes reproduce diploid sporophytes by special reproductive cells (e.g., mitospores).I note, however, that it is also known that "haploid sporophytes" reproduced by haploid gametophytes and "diploid gametophytes" reproduced by diploid sporophytes are observed in empirical conditions [58,59].

Population dynamics without density dependence
First, the population dynamics of a haploid-diploid population without density dependence is analyzed by updating the matrix model [28].The haploid (H) and diploid (D) densities are forecasted to time t+1 from their realized densities at time t and their vital rates: Here, the survivorship of haploid and diploid reproductive cells (tetraspores and carpospores for red algal species) while suspended in the water column is γ H and γ D , respectively; the fecundity of the haploid and diploid stages per individual is b H and b D , respectively; and the mortality rates of haploids and diploids are m H and m D (adult survival is 1-m), respectively.The b H is the fecundity rate of haploid females (production of carposporophytes for red algal species).However, not all haploids are females as some are males.Assuming a 1:1 sex ratio, 1/2 of the haploids are females (or we can also assume the species is monoecious).Hence, in the matrix model, haploid fecundity is b H /2.Then, comes the reproductive cost (the probability of fertilization success) f, leading to fb H /2. In the following analysis, I assume this cost (cost of sex) is defined as σ = f/2.Fractions a H and a D of the reproductive outputs of the haploid and diploid stages are asexual, skipping the sexual loop and developing directly into haploids and diploids (sporogenesis).When a H = a D = 0, Eq (1) describes the case of a fully sexual (obligate sexual) system [28].The dominance patterns of two stages following spatial or seasonal patterns are often observed in the field, however for tractability, constant life history parameters and a well-mixed population are assumed.

Population dynamics with density dependence
A density-dependent population of macroalgae that consists of haploids and diploids can be described by the following differential equations:

Ploidy ratios with a density-independent model
With a density-independent model, the population size increases exponentially.In this model, the right leading eigenvector, which gives the stable distribution of haploid and diploid abundances in the population, is conceptually equivalent to the stable age distribution in the Leslie matrix model.The eigenvalues and right eigenvectors for the full model are provided in the S1 File.The haploid frequency in the population, ρ H , can be calculated as follows: ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi where the fitness components of haploids and diploids are defined as, and the superscripts S and A refer to sexual and asexual reproduction, respectively.Note that the fitness components for sexual reproduction (w S H and w S D ) is defined by the mortality of different ploidy individuals.This is because the reproductive cells for sexual reproduction reproduced by individuals of focal generation (e.g., haploid) will develop to the opposite generation (e.g., diploid), so mortality should be different from the focal individual.
If the species has an "ideal" isomorphic life cycle and exhibits full sexuality (obligate sexuality), the number of parameters can be reduced as follows: and a D = 0.In this case, the haploid ( � H) to diploid ( � D) ratio in asymptotic density is, If the cost of fertilization is ignored (f = 1), then σ = 1/2 and Eq (4) is the same as the classical ploidy ratio of ffi ffi ffi 2 p : 1 ([28], p. 1664).If the fecundity of haploids is assumed to be twice that of diploids, because the haploid cell size tends to be half the diploid cell size, then γ = γ H = γ D , m = m H = m D , and b H = 2b D , and the ploidy ratio is, ffi ffi ffi ffi ffiffi 2s If the cost of fertilization is ignored (f = 1), then σ = 1/2 and the ploidy ratio in this case is 1:1.

Ploidy ratio of a fully sexual species with density dependence
In the density-independent model, population size increases infinitely, which is unrealistic.To avoid this problem, the density-dependent model is analyzed next.For tractability, the model for a species that exhibits full sexuality (a H = a D = 0) is analyzed first.In this case, strict solutions for both biomass (T ¼ Ĥ þ D) and haploid frequency (r H ¼ Ĥ=ð Ĥ þ DÞ) are obtained at population equilibrium ( Ĥ; D), Here, the fitness of haploids and diploids is defined as w H = (σb H γ D )/m D and w D = (b D γ H )/ m H , respectively.Note that the mortality rate in the fitness of haploids (w H ) should be defined by the diploid mortality (m D ), and vice versa.This is for exactly the same reason as the definition of the fitness components in the full model (w S H and w S D ).This haploid frequency (ρ H in Eq 6) is equivalent to the haploid-diploid equilibrium frequency in the Wright-Fisher model in Bessho and Otto ([60]; their Eq (D2)).Furthermore, the haploid frequency, ρ H , can also be represented as follows: r H ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ðb D g H Þ=m H p ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ½ðb D g H Þm D Þ p ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi This haploid frequency is the same as that in the Moran model in Bessho and Otto ([60]; their Eq (A8b)).
Note that when the mortalities of haploids and diploids are equal, r H =r D ¼ ffi ffi ffi ffi ffi ffi , which is the gametophyte-to-sporophyte ratio in Fierst et al. ( [30], their Eq (15))

Parameter dependence in species with obligate sexuality
The basic reproductive number of species with obligate sexuality is the geometric mean of haploid fitness and diploid fitness, R 0 ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi w H w D p (Appendix C2 in S1 Appendix).Thus, the equilibrium biomass calculated using this measure is, Thus, the biomass (total population size) holds constant at a positive value at equilibrium when R 0 >1.Furthermore, when R 0 >1, the trivial equilibrium, (0, 0), is unstable and the haploid-diploid population exists stably as in Eq (6) (Appendixes B3 and B4 in S1 Appendix).
Fig 2 illustrates the parameter dependence of the haploid frequency ρ H and equilibrium biomass, scaling the population size by δ in Eq (7) (i.e., T ¼ Td ¼ R 0 À 1; hereafter, the biomass measure).The ploidy ratio and biomass both depend on haploid fitness (w H ) and diploid fitness (w D ).Because haploids reproduce diploid individuals and vice versa, haploids are dominant in a population with higher diploid fitness.Note, however, that fitnesses are defined by the mortality of opposite ploidy individuals (see the definition).For example, a high haploid mortality rate (high m H ) reduces the diploid fitness (low w D ) and consequently increases the frequency of diploids (low ρ H ).

Ploidy ratio in asexually reproducing species under density dependence
For asexually reproducing species, strict solutions for biomass (T ¼ Ĥ þ D) and haploid frequency (r H ¼ Ĥ=ð Ĥ þ DÞ) are obtained at the population equilibrium ( Ĥ; D) (see details in the S1 File), ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi https://doi.org/10.1371/journal.pone.0295409.g002 ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi where I define the fitness components of haploids and diploids as follows: When the fitness components of asexual reproduction are equal between haploids and diploids, w A H ¼ w A D , haploid frequency is simply, Interestingly, the ploidy ratio in the density-dependent model, Eq (8B), is equivalent to that in the density-independent model (3) when the mortalities of the haploid and diploid stages are equal (m H = m D ).
If equal survivorships of the reproductive cells are assumed (γ = γ H = γ D ) along with equal mortality rates (m = m H = m D ), then the reproductive parameters can be simplified and the ploidy ratio at equilibrium becomes ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi where bH ¼ ðb H gÞ=m and bD ¼ ðb D gÞ=m.Eq (10) is the equilibrium ploidy ratio in Bessho and Otto ([53], their Eq (A.4)).Given that they were using the Wright-Fisher model and thus ignored differences in mortality between haploids and diploids, this match makes sense.

Parameter dependence of asexually reproducing species
For tractability, consider first the case where the fitness components of haploids and diploids are symmetrical for both sexual and asexual reproduction, w S ¼ w S H ¼ w S D and w A ¼ w A H ¼ w A D .Under this condition, the basic reproductive number is the sum of the fitness components of sexual and asexual reproduction, R 0 = w A +w S .When this measure is larger than one, the trivial equilibrium is unstable, and the nontrivial equilibrium is stable (Appendix C1 in S1 Appendix).
If the assumption of symmetry is relaxed, the basic reproductive number (R 0 ) becomes somewhat more complicated (see Eq (C.6)), but Eq (7) still hold.Unfortunately, it is not possible to simplify the stability condition, so it was checked numerically.It can be concluded that this haploid-diploid system never shows bistability and oscillation.When the basic reproductive number is larger than one (R 0 >1), the trivial equilibrium becomes unstable and the dynamics approaches a nontrivial equilibrium in which the haploid frequency is described by Eq (8B).
Figs 3-7 illustrate the parameter dependence of the haploid frequency (ρ H ) and the biomass measure ( T ¼ R 0 À 1).Both the ploidy ratio and biomass depend on the fitness components (w S H ; w S D ; w A H , and w A D ).Because asymmetry between haploids and diploids (i.e., the cost of sex, σ) is absorbed into the transformed parameter w S H , the relationship between haploids and diploids becomes symmetrical when comparing the w's.
The parameter dependence of the fitness components for sexual reproduction (Figs 3 and  4) can be understood by comparing to the case of obligate sexuality (Fig 2).As shown by Eq (9), when the asexual fitness of the species is symmetrical between haploids and diploids (w A H ¼ w A D ), the parameter dependence is similar to that for an obligate sexual species (Figs 3B, 3E, 4B, and 4E).When the asexuality of haploids is stronger than that of diploids (w A H > w A D ), then the region in which the haploid fraction is dominant in the population (ρ H >0.5) becomes larger.The reverse (w A H < w A D ) is also true.Interestingly, the haploid frequency in a population, ρ H , depends on the difference in the fitness components for asexual reproduction, w A H À w A D , when w S H and w S D are fixed.Hence, the slopes of the haploid frequency curves become one (Fig 5A -5C).This result means that the ploidy ratio in a haploid-diploid population can be described by three parameters (see Eq (8B), which depends on the sexual fitnesses only through w A H À w A D ).When the haploid frequency is illustrated as a function of the difference in the fitness components for asexual reproduction (w A H À w A D ), the curves are sigmoidal (Fig 5G).When the fitness components for asexual reproduction of haploids and diploids are similar (w A D � w A H ), then the basic reproductive number can be approximated by the sum of the geometric mean of the fitness components for sexual reproduction and the arithmetic mean of the fitness components for asexuality, R 0 � ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi )).Hence, when the fitness components for asexual reproduction are given, the biomass measure ( T ¼ R 0 À 1) is a function of the geometric mean of the fitness components for sexual reproduction (Fig 4).

Ploidy ratio in a haploid-diploid population
Mixtures of free-living haploid gametophytes and diploid sporophytes are widely observed in macroalgal populations.This article presents a general haploid-diploid model that explicitly accounts for density dependence and asexual reproduction and can be used to obtain an exact solution for the ratio of haploid gametophytes to diploid sporophytes in a population at steady state.The analysis (Eq 8B) reveals that the proportions of gametophytes and sporophytes in a wild population are determined by the balance of four fitness components: sexual reproduction of haploids (w S H ), sexual reproduction of diploids (w S D ), asexual reproduction of haploids (w A H ), and asexual reproduction of diploids (w A D ).Furthermore, the difference in the asexual demographic force (w A H À w A D ) is important.Because all life history parameters (i.e., fecundity, mortality, survivorship of reproductive cells, and the cost of sexuality) are embedded in these fitness components, the composition of a haploid-diploid population is entirely determined by three parameters (w S H ; w S D , and w A H À w A D ).I note that the ploidy ratio that is obtained is robust to changes in the mode of competition (Appendix D in S1 Appendix).
The result is consistent with the ploidy ratios reported by previous studies for the special cases considered.For example, the classical haploid-to-diploid ratio of Thornber and Gaines [28], ffi ffi ffi 2 p : 1, is consistent with Eq (8B) when the life history parameters of two stages are assumed to be identical ("ideal" isomorphy).The equilibrium ploidy ratios in the Wright-Fisher and Moran population genetic models [53,60] also correspond to a special case of Eq (8B).
In summary, the model presented here generalizes the findings of previous studies [28,30,53,60] and mathematically demonstrates the intuitive conclusion that the frequency of haploidy in a population is determined by the balance of fitness components, with sexual reproduction in a haploid or diploid individual leading to an increase in the other's frequency and asexual reproduction in a haploid or diploid leading to an increase in its own frequency.It should be noted, however, that the fitness component for sexual reproduction in our analysis is defined by the mortality rates of the opposite ploidy stage.For example, an increase in haploid mortality decreases the fitness component of haploid for the asexual reproduction while simultaneously decreasing the fitness component of diploid for the sexual reproduction.

Stability and biomass of haploid-diploid populations
Despite the fact that two-dimensional systems of differential equations with density dependence may exhibit a variety of dynamics, including bistability and oscillations, previous research has lacked studies of population stability (e.g., [61]).With the help of numerical calculations, the dynamics of a haploid-diploid population has been found to exhibit a stable equilibrium state without oscillations in the parameter region in which the population exists (positive biomass).
Furthermore, the persistance of a population and the equilibirum population size depend critically on the basic reproductive number, R 0 (Eq (C.6)).In particular, when the fitness components for asexual reproduction of haploids and diploids are similar (w A H � w A D ), then Thus, the basic reproductive number is approximately the sum of the geometric mean of the fitness components of haploids and diploids for sexual reproduction and the arithmetic mean of the fitness components of haploids and diploids for asexual reproduction.Hence, the contribution of sexual and asexual reproduction to the existence of a population can be estimated by calculating these geometric and arithmetic means of the fitness components.
Empirical studies have shown that some macroalgal populations are maintained by asexual reproduction (e.g., [51]).Interestingly, the model presented here suggests that, when the sexual reproduction term is smaller and can be ignored, the basic reproductive number can be approximated as ( These results suggest that when asexual reproduction is dominant in a population, the ploidy with the larger fitness component mainly determines population persistence and total population size.

Ecology and evolution of haploid-diploid life cycles
The life cycle diversity observed in macroalgae has been of interest in terms of the evolution of haploid versus diploid [9][10][11]13].However, a simple population genetic model proposed the problem that life cycles in which both stages are stably maintained (biphasic life cycle) do not evolve, and the importance of ecological effects began to attract attention [57,61,62].As an influential study, Hughes and Otto [61] found that ecological differences between two ploidy stages can lead to the evolution and maintenance of biphasic life cycles.In contrast, the evolution of variety in haploid-diploids (i.e., heteromphy versus isomorphy, and haploid-dominant versus diploid-dominant) is still obscure.Because the evolution in a haploid-diploid should be analyzed the mutant dynamics in a resident population that reached the equilibrium (e.g., [53]), the stable ploidy ratio in this article contributes to this analysis as the first step.
The haploid-diploidy in macroalgae has raised another question for researchers: how is the ploidy ratio in the field determined?The relative importance of the difference in life history parameters (e.g., fertility, survivorship, and fertilization success) between haploids and diploids was analyzed by simple [28,30] or complex [31,38] mathematical models.However, our study reveals that essentially life history parameters can be summarized in only three components; sexual reproduction of haploids (w S H ), sexual reproduction of diploids (w S D ), the difference in the asexual demographic force (w A H À w A D ).Of course, the measurement of individual life history parameters remains important, but in the future, when classifying and understanding the diversity of ploidy ratios observed in the field, the question of how the three components are balanced is more essential than the question of which life history traits are important.Furthermore, if the difference in asexual demographic forces is negligible, the balance between the geometric mean of sexual demographic forces ( ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi w S H w S D p ) and the arithmetic mean of asexual demographic forces (ðw A H þ w A D Þ=2) will also be an important property in applications, in the sense that population stability and biomass can be evaluated by the sum these two means.

Future perspectives
This article presents an analysis of the ploidy ratio and stability of macroalgal populations observed in the field.However, several issues remain to be addressed.First, it is assumed in this paper that density dependence is regulated by the population biomass.However, asymmetry between gametophytes and sporophytes may also affect the mode of competition.In other words, the question is whether the introduction of ploidy-specific competition coefficients (e.g., [61]) would change the conclusion.Similarly, since this model assumes that the parameters that determine the strength of the density dependence (δ) are equivalent between haploids and diploids, the difference of δ between ploidy phases is also an interesting problem.Because R 0 is a measure of spread when rare determining (stability of the trivial equilibrium), it is not affected by the nature of density dependence.However, it remains unclear how the stability of a nontrivial equilibrium is affected by the competition coefficient and the strength of the density dependence; thus, it is important to extend the model to investigate these points.
Second, although it is assumed that all life-history parameters are constant, wild populations often exhibit seasonal oscillations.This issue is also related to the stability of population dynamics in the first point.It is an important task to determine whether oscillations are an intrinsic property of the population dynamics or whether they are due to parameter changes caused by seasonality and to investigate the effects of parameter changes on population dynamics.
The demography and the ploidy structure of a population are fundamental to understanding the evolution of diverse and complex life cycles.The analyses and formulas developed here provide a basis for analyzing ecological data to discuss demography and evolution in haploiddiploid species.

Fig 1 .
Fig 1. Diagrams of macroalgal life cycles.(a) A diagram of the reproduction process assumed in the model.(b) A diagram of the alternation of generations between haploid gametophyte and diploid sporophyte (haploid-diploid life cycle).In this diagram, monoecious is assumed.(c) A diagram of the general life cycle of Florideophyceae (triphasic life cycle).(d) A diagram of the asexual reproduction in macroalgae.Red indicates the haploid phase and blue indicates the diploid phase.The circled states are free-living stages.https://doi.org/10.1371/journal.pone.0295409.g001

Fig 2 .
Fig 2. The haploid frequency ((ρ H ¼ Ĥ =ð Ĥ þ DÞ) and the biomass measure ( T ¼ R 0 À 1) in a population of a species with obligate sexuality (a H = a D = 0).(a) Haploid frequency curves; the red line represents the case that the haploid fraction constitutes exactly half of the population (ρ H = 0.5).(b) Curves for the biomass measure.The red line represents the case that R 0 is exactly equal to one ( T ¼ 0).Population extinction occurs in the region where T is negative.(c) The haploid frequency as a function of haploid fitness (w H ) for different values of diploid fitness (w D ).Dashed lines indicate values for extinct populations ( T < 0).(d) The biomass measure as a function of haploid fitness (w H ) for different values of diploid fitness (w D ).

Fig 3 .
Fig 3. Haploid frequency (ρ H ¼ Ĥ =ð Ĥ þ DÞ) in a population of a species with asexual reproduction when the fitness components for asexual reproduction (w A H and w A D ) are given.(a),(b), and (c) Haploid frequency curves.The red lines represent the case where the haploid fraction constitutes exactly half of the population (ρ H = 0.5).(d),(e), and (f) The haploid frequency as a function of the haploid fitness component for sexual reproduction (w S H ) for different values for the diploid fitness component for sexual reproduction (w S D ).Dashed lines indicate values for extinct populations ( T < 0).The given parameter values are (a),(d) w A H ¼ 0:8 and w A D ¼ 0:2; (b),(e) w A H ¼ 0:5 and w A D ¼ 0:5; and (c),(f) w A H ¼ 0:2 and w A D ¼ 0:8.https://doi.org/10.1371/journal.pone.0295409.g003

Fig 4 .Fig 5 .
Fig 4. The biomass measure ( T ¼ R 0 À 1) in a population of a species with asexual reproduction when the fitness components for asexual reproduction (w A H and w A D ) are given.(a),(b), and (c) Curves for the biomass measure.The red lines represent the case that R 0 is exactly equal to one ( T ¼ 0).(d), (e), and (f) The biomass measure as a function of the haploid fitness component for sexual reproduction (w S H ) for different values of the diploid fitness component for sexual reproduction (w S D ).The given parameter values are the same as in Fig 3. https://doi.org/10.1371/journal.pone.0295409.g004

Fig 6 .
Fig 6.The biomass measure ( T ¼ R 0 À 1) in a population of a species with asexual reproduction when the fitness components for sexual reproduction (w S H and w S D ) are given.(a),(b), and (c) Curves for the biomass measure.The red lines represent the case that R 0 is exactly equal to one ( T ¼ 0).(d),(e), and (f) The biomass measure as a function of the haploid fitness component for asexual reproduction (w A H ) for different values of the diploid fitness component for asexual reproduction (w A D ).The given parameter values are the same as in Fig 5. https://doi.org/10.1371/journal.pone.0295409.g006

Fig 7 .
Fig 7. The biomass measure ( T ¼ R 0 À 1) in a population of a species with asexual reproduction when the fitness components for sexual reproduction (w S H and w S D ) are given.(a),(b) Curves of the biomass measure.(c),(d) The biomass measure as a function of the haploid fitness component for asexual reproduction (w A H ) for different values of the diploid fitness component for asexual reproduction (w A D ).The given parameter values are (a),(c) w S H ¼ 0:5 and w S D ¼ 0:5, and (b),(d) w S H ¼ 3:0 and w S D ¼ 3:0.https://doi.org/10.1371/journal.pone.0295409.g007 Here, the mortality rate is assumed to increase with the biomass of the population (H(t)+D(t)) at rate δ.When population dynamics are modeled with density dependence, other modes of action, besides increasing mortality, should be considered (e.g., decreasing fecundity).Other types of density dependence are considered in Appendix D in S1 Appendix.